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ABSTRACT 

The aim of this paper is to demonstrate the accuracy of our knowledge of close stellar passage 
distances in the pre-GAIA era. We used the most precise astrometric and kinematic data avail¬ 
able at the moment and prepared a list of 40 stars nominally passing (in the past or future) 
closer than 2 pc from the Sun. We used a full gravitational potential of the Galaxy to calculate 
the motion of the Sun and a star from their current positions to the proximity epoch. For this 
calculations we used a numerical integration in rectangular, Galactocentric coordinates. We 
showed that in many cases the numerical integration of the star motion gives significantly dif¬ 
ferent results than popular rectilinear approximation. We found several new stellar candidates 
for close visitors in past or in future. 

We used a covariance matrices of the astrometric data for each star to estimate the accur¬ 
acy of the obtained proximity distance and epoch. To this aim we used a Monte Carlo method, 
replaced each star with 10000 of its clones and studied the distribution of their individual 
close passages near the Sun. We showed that for contemporary close neighbours the precision 
is quite good but for more distant stars it strongly depends on the quality of astrometric and 
kinematic data. Several examples are discussed in detail, among them the case of HIP 14473. 
For this star we obtained the nominal proximity distance as small as 0.22 pc 3.78 Myr ago. 
However there exist strong need for more precise astrometry of this star since the proximity 
point uncertainty is unacceptably large. 
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1 INTRODUCTION 


Soon after the Oor dl950l) published a paper on the existence 
of the comet cloud surrounding the Sun, many authors started 
to investigate the influence of real nearby stars on cometary or¬ 
bits. One of the earliest complete (at that mo ment) search fo r 
stellar passages near the Sun was published by iMakox^ l ll964h . 
He listed all past stellar proximity epochs and d istances based 
on the first edition of the Gliese Catalogue (Il957h and obtained 
with a rectilinear approximation. In recent twenty five years sev¬ 
eral papers were published in this field, among others Matthe^ 
I 994 I; Mullari & OrlovI 199^ Dvbcz vnski & Kankiewi^ 19991; 


Garcia-Sa nchez et all I1999L 120011; iDvbczvhski 200d ; Bobvle^ 

2010hl lal; Ijimenez-Torres et al.ll201 ll . The main reason for a con¬ 
tinuous interest in determining close star passages is the pro¬ 
gress in obtaining stellar data, especially proper motions, paral¬ 
laxes and radial velocities. Additionally, investigations of long 
period comets past and future dynamics are more and m o re de¬ 
tailed (see for example Krolikowska & Dvbczvnskil 201^ 20G 


Dybczyhski & Kr6 1ikowskall201 ll . l2015l ; iKrolikowska et ^ 2012 ; 

Krolikowskal 2014 ). As a result best possible knowledge on 
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S iotential stellar pertu rbers of such motion is of great interest 
Fouchard et alj2oTII) . 

With the publication of the first Hipparcos catalogue i ll9971 . 
hereafter HIPl) the number of stars with known para llaxes has in¬ 
creased significantly. Hipparcos mission ( IESA||19^ was a mile¬ 
stone in our last decades advance in gathering data on the spatial 
distribution of stars in the solar neighbourhood. This mission gave 
us precise positions, parallaxes and proper motions for -120 thou¬ 
sands of stars from carefully and i n advance selected list, form¬ 
ing the Hipparcos Input Catalogue jTuron et alJI 19921) . Hipparcos 
catalogue is c omplete up to ~ 25 pc for all stars brighter than 
My = 9 mag dlahreiB & Wielenlll997l) and for the massive stars 
in this population the completeness reaches a much larger distance. 
There are two main reasons for which it is difficult to obtain pre¬ 
cise information on the spatial velocity for all Hipparcos stars. First 
- the Hipparcos proper motions, while of great formal precision, 
are based on relatively short time interval, therefore for a large 
number of stars the y differ from the mean, secular proper motions 
dWielen et alJI 19991) . Several attempts were made to combine Hip¬ 
parcos proper motions with those from a ground based observa¬ 
tions, see for example dWielen et alJll999| ; iHoogerwerf & Blaauwl 
l200d;IWielen et alj200lF The most fruit-full attempt was the con¬ 
struction of Tycho-2 catalogue dH0g et alj|2000al f9) however com¬ 
bining Tycho-2 proper motions with parallaxes from HIP-1 intro¬ 
duces some inconsistency in treating astrometric data. Recently a 
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new reduction of Hipparcos raw data were performed and a second, 
significantly improved ver sion of the Hipparcos catalogue (HIP2, 
Ivan Leeuwenll20075 . 1201 ll ) have been published. It improves sig¬ 
nificantly precision of both proper motions and parallaxes. 

As it concer ns radial velocity, se veral large projects 
(see for examp le Grenier et al.l ll999al S iNidever et al.l l2002l : 
iNordstrbm et alj|2004tl increased the number available measure¬ 
ments and a large Pulkovo compilatio n of radial v e locitie s of 35 
495 Hipparcos stars was published by iGontcharovI ( l2006h . Using 
HIP2 catalogue, Pulkovo compilation of radi al velocities and sev¬ 
eral ot her sources of astrophysical parameters I Anderson & Francis! 
(l2012h published an extended Hipparcos compilation of stellar 
data, known as the XHIP catalogue. 

The second reason for investigating stellar passages is our 
increasing ability to calculate the stellar path relative to the Sun 
with increasing accuracy, taking into account the Galactic grav¬ 
ity field. Using such advanced methods we should additionally ask 
about the accuracy of the results and their dependency on the used 
data. Authors of several previous papers also estimated the accur¬ 
acy o f their results in more or les s approximate manner. For ex¬ 
ample iGarcfa-Sanchez et al.l ( 1200 ih simply used the root of a sum 
of squared error contributions from two components of proper mo¬ 
tions, parallaxes and radial velocities. The Monte Carlo method of 
estimating the accuracy of the distance of a stell a r clos e passage 
near the Sun was first used by iMullari & OrlovI d 19961) but in a 
simplified manner (for example with assumed 3 kms^* error for 
all radi al velocities). Simi lar attempt was perfo r med b y fBobvlevI 
(I2010bl l3) and recently by Ijimenez-Torres et ^ d 201 ll) but in all 
cases only formal errors were taken into account, ignoring mutual 
correlations between astrometric parameters. 

In the present work we used astrometric parameters taken dir¬ 
ectly from the HIP2 catalogue instead of XHIP mainly because 
HIP2 presents also a covariance matrix which is necessary to apply 
the advanced method for the accuracy assessment. We used radial 
velocities from the XHIP catalogue. For several selected stars we 
presented comparison of the results based on astrometry and radial 
velocities from some other sources. 

This paper provides the most up to date information on the 
closest stellar approaches to the Sun, obtained from latest astromet¬ 
ric data and radial velocities augmented with the more elaborated 


Monte Carlo assessment of the accuracy of the obtained minimal 
distances and their epochs, based on full covariance matrices in¬ 
cluded in the HIP2 catalogue. Of course we expect significant im¬ 
provements from the Gaia mission therefore our computer codes 
are fully prepared to include new data. 

In the next section we describe methods of our calcula¬ 
tions. Section [3] consists of the results obtained for selected stars. 
In section we prese nt a discussion with a very recent paper 
by ISailer-JonesI ( 1201 4 hereafter B-I), indpendently prepared and 
made available to the public when this work was finished. In the 
last section some conclusions and prospects are drawn. 


denotes the multiplication by cos 5 ) and radial velocity ly. The 
starting data for a numerical integration, namely position and velo¬ 
city components in the Galactocentric frame are to be obtained as 
follows. We first calculate the heliocentric distance of a star from 
the formula: 


ni = 


1000 

It 


( 1 ) 


where the parallax 7t is expressed in miliarcseconds [mas] 
wh at gives ri, in parse cs. According to comments expressed 
by Ivan LeeuwenI ( l2007al . chapter 3, page 86 ) we resisted here 
to add the so called Lutz -Kelker bias correction, proposed by 
lAnderson & Francis ( |2012|) . 

Next, heliocentric, equatorial rectangular coordinates of a star 
(in parsecs) are given by: 


Xh = r/, ■ cos a cos d 

yh = fh ■ sin a cos d ( 2 ) 

Zh = 0 , ■ sin 5 

If proper motions are given in mas/yr and radial velocity is 
expressed in kms^* than the heliocentric, equatorial rectangular 
velocity components (in pc/Myr) might be calculated as follows: 


VI = s- ri,- flat/cos5 
V2=s-rh-lls 

V 3 =k-Vr 

X;, = — COS 5 sin a • vi — sin 5 cos a-V2+ cos d cos a ■ V 3 
yij = cos d cos a • V] — sin 5 sin a • V 2 -F cos 5 sin a • V 3 
th = cos 5 • V2 -F sin 5 • V3 

where s = 0.0048481368 is a coefficient for angular units con¬ 
version and k = \ .022689369 is a coefficient for velocity units con¬ 
version. Throughout this paper we use the following units: parsec 
[pc] as the distance unit, solar mass [Mq] as the mass unit and 
million of years [Myr] as the time unit. The constant of gravity ex¬ 
pressed in these units is: G= 4.498297316 x 10^^. 

Reorientation of the equatorial frame into Galactic one (still 
heliocentric) involves three rotations: 

rh.gai = Rz{90° - 0)RA9O° - 5„)R^(90° -F ao)rh.equ 

where 9 = 122.°93191857, a„ = 12''51'"26?27549 and do = 
27°07'42"7043 are the positional angle and north Galactic pole 
equatorial coordinates defining the Galactic frame orientation while 
Rj and Rj denot e the rota t ion w ith respect to OX and OZ axes 
respectively. See lUiu et al.l ( 1201 ih for a recent discussion on the 
Galactic frame orientation. 

To move the origin from the Sun to the Galactic centre it is 
necessary to add the Galactic position and velocity of the Sun: 


2 METHODS OF CALCULATIONS 
2.1 Units, definitions, reference frames 

To achieve the aim of this work we have to study stars (including 
the Sun) motion under the gravitational influence of the Galaxy. 
This can be performed only for stars with the full 6 D data available, 
typically expressed in the equatorial frame, i.e. right ascension a, 
declination S, parallax 7t, proper motions fig (a star subscript 


R — ^h.gal ~F R0, R — ^h.gal ~F Rq , 

The values for the solar position Rq = (xQ,yQ,Z0) and ve¬ 
locity R0 = (m,v,w) components should be chosen in accordance 
with the adopted Galactic potential model. In our calculation we 
use R 0 = (j^ 0 ,y 0 ,Z 0 ) = (—8400,0,17) in pc. Since the vertical 
position of the Sun with respect to the Gala ctic disk plan e is still un¬ 
certain we decided to follow arguments of lloshil ( l2007h and adopt 
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zq =+17 pc. The reader should be warned that during this transla¬ 
tion we kept the orientation of all axes, so in a Galactocentric frame 
the OX axis is directed opposite to the Sun direction, the OY is dir¬ 
ected in accordance to the Sun rotation around the Galactic centre 
and the OZ axis is directed to the north Galactic pole, what makes 
this system a right-handed one. 

As i t concerns the Sun p eculiar velocity, we use values pro¬ 
posed bv lirrgane et alj ( I2Q13I) . what after adding the Sun Galactic 
rotation velocity resulting from the adopted potential model gives: 

Rq = {u,v,w) = (d-l 1.1,+254.24,-f7.25) 

= (+11.352,+260.011,+7.41) pc/Myr. 

2.2 Drawing a stellar clone 

In order to estimate the uncertainty of the close stellar passage para¬ 
meters we replace each considered star with a large number (typ¬ 
ically 10 000) of its clones, drawn from a multivaria te normal dis¬ 
tribution. Thi s is possible because in HIP2 catalogue lvan LeeuwenI 
(l2007bLl201 ih included full information on the covariance matrix of 
astrometric parameters. Such a procedure possesses an evident su¬ 
periority over the individual, independent ran dom draw i ng of ea ch 
parameter (which was use d for example by ISobvlevI ( I2010al f9): 
Ijimenez-Torres et aP (1201 ih ) what ignores their obvious mutual de¬ 
pendence. We demonstrate this at the end of this section. There are 
no radial velocities in HIP2 so we have to draw radial velocity inde¬ 
pendently using data from XHIP catalogue. Unfortunately for some 
stars the formal error is unknown and the authors insert 999 km s^* 
in catalogue. One of such stars is HIP 21539 for which we adopted 
30 kms^* as formal e rror. We used this value bec ause this radial 
velocity measurement dSarbier-Brossat et al.lll994) most probably 
comes from observation with an obiective prism which have the 
precision of this order. 

To generate random vectors x from a multivariate normal dis¬ 
tribution (m, £) with a given vector of means m an a given co- 
variance matrix Z one can proceed as follows: 

(i) Decompose the given covariance matrix (e.g. with Cholesky 
method) to obtain the matrix G such, that Z = GG^. 

(ii) Compose a vector y of n independent random scalars y, from 
a standard normal distribution: y,- ~ ^(0,1 ). 

(iii) Obtain vector x from: x = m + Gy. 

Here are some details on applying the above recipe for a HIP2 star 
cloning. In what follows we assume the 5 parameter solution but 
the analogous approach can be applied to stars with the 7- or 9- 
parameters solutions. For each star in HIP2 catalogue given is an 
upper triangular ’weight matrix’ U. According to the catalogue de¬ 
scription it should be related to the covariance matrix Z by a simple 
relation: 

Z^*=U^U (3) 

However this is not exactly true (van Leeuwen, 2013, personal 
communication). To obtain matrix U from catalogue data it is ne¬ 
cessary to re-scale its diagonal elements (mj ,M 3,M6 ,mio and ni5 from 
catalogue) multiplying each of them by a factor of (c7o/Ax,) where 
Axi is one of Aa*, A5, A;r, Apa* and Apg, i.e. given in the cata¬ 
logue formal errors of the astrometric parameters and Oq is a ref¬ 
erence variance (in other words: a ’unit weight observation uncer¬ 
tainty’) for a given star to be cloned. To calculate Oq value one 


should retrieve three additional numbers from the stars’ record in 
the catalogue: ’goodness of fit’ - F2, ’number of field transits used’ 
- Ntr and ’percentage of rejected data’ - FI and use the following 
relations: 


dof =yVfr(l-Fl/100)-5 

(4) 

c =2/{9-dof) 


X =(l-c + F2-f^fl^-dof 

(5) 

O’o =X/'/ddf 



The above is a reverse of the procedure used for a ’goodness 
of fit’ calculation (van Leeuwen, 2013, personal communication). 
When calculating the ’degrees of freedom’ (dof, equation|4} again 
the 5-parameters solution is assumed here. It is worth to note that 
formal errors of all five astrometric parameters, namely: Act*, AS, 
Alt, Apa* and Apg presented in the catalogue are calculated (van 
Leeuwen, 2013, personal communication) by multiplying corres¬ 
ponding diagonal elements from the covariance matrix by the ref¬ 
erence variance Oq . 

Summarising all the above we obtain the matrix U in a form: 


obiii/Aa* 

U2 

«4 

u-i 

«11 

0 

fToU3 /AS 


Wg 

»12 

0 

0 

aoUi/An 

Ug 

Hl3 

0 

0 

0 

Ob»lo/Apa* 

1(14 

0 

0 

0 

0 

ob«i5/Apg 

(f 


According to equation[^we have: 


Z= (U^U)^' =U^'(U^)^' = GG^ (7) 

what shows, that G = is suitable for our purpose of draw¬ 
ing a star clone according to the recipe presented at the beginning 
of this section. 

We tested how different methods of drawing clones affect the 
dispersion of the proximity position. To do this we repeated our cal¬ 
culation for all selected stars by drawing all parameters of stellar 
clone independently. In most cases distribution was only slightly 
different in shape, but in few cases, for example HIP 25240 (see 
Fig.B dispersion was significantly higher than when we used co- 
variance matrix. This shows that the use of the covariance matrix, 
apart from being closer to the data, in some cases can significantly 
improve the accuracy of the close stellar approach determination. 

2.3 Galactic motion of stars 

To study the Galactocentric motion of a star we use a numerical 
integration of its equations of motion expressed in rectangular co- 
or dinates, utilising the well known, fast and accurate RA15 routine 
bv lEverhartI ( Il985l) . To describe the Galac tic gravitational potential 
we use Model I from lirrgang et al.l l l2013l) . Since we are interested 
in motion of nearby stars from their current positions backward or 
forward to the moment of their closest heliocentric position we deal 
with rather small time intervals of order of 10 Myr and small helio¬ 
centric distances, not exceeding 200 pc, with the only one exception 
of the star HIP 33369 with the current heliocentric distance of 424 
pc (such a large distance makes our results for this star completely 
unreliable as it is shown in Section[3](- For such a small solar vicin¬ 
ity a more sophisticated Galactic potential model, accounting for 
non-spherical central bulge or for spiral arms is not necessary, as 
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Figure 1. Differences in a dispersion of clones of HIP 25240 sampling from six one-dimensional Gaussians (red dots) and drawing clones from covaiiance 
matrix (black dots). All clones are stopped at the closest proximity epoch. Each green point mark the nominal result. Blue circle represents here the boundary 
of the Oort cloud at 0.5 pc. 


Table 1. Model I parameters from lirrgang et al.l i2013h 


Parameter 

Value 

the distance of the Sun from the Galactic centre /?0 

8400 pc 

Galactic bulge mass 

9.51x1O‘’M0 

Galactic disk mass 

66.4x IO’Mq 

Galactic halo mass M/, 

IS.VxIO’Mq 

bulge characteristic distance bf, 

230 pc 

disk characteristic distance a^j 

4220 pc 

disk characteristic distance 

292 pc 

halo characteristic distance «/, 

2562 pc 

Galactic halo cut-off parameterA 

200000 pc 

Galactic halo exponent parametery 

2 (fixed) 

Galactic disk matter density near the Sun po 

0.102 Mq/pc^ 

Galactic rotational velocity of the LSR Vg 

242kms-‘ 


^h{R) = ' 


^ I TT^-rln 


Oh I r-1 


x(r-i) 


\ (r-i) 






x(r-i) 


1 + 


Mh \‘‘hj 


R 




\(r-i) 


ifR<A 


elsewhere. 


what, after adopting 7=2 and choosing the first equation (we 
certainly do not go as far && R = A =200 kpc!) reduces to: 


’^h{x,y,z) 


^ A ( fl/i ++ 
ail ( ( fl;,+A 


ah+Aj 


The equation of motion of a single point mass (a star) under 
the potential described above, expressed in a rectangular Galactic 
(and Galactocentric) frame are; 


was checke d for example by Garcia-Sanchez et"^ (1200 ih or more 
recently bv Ijimenez-Torres et alj ( 1201 ih . 

Gravitational potential 'h(r,z) considered here is the sum 
of a central bulge component ^h{R) (spherically symmetric), an 
axisymmetric disk <f>^(r,z) and a massive spherical Galactic halo 
^h{R) (dark matter included): 


where {r,(p,z) are Galactocentric cylindrical coordinates and 
R = ^yp+z^ is a Galactocentric spherical radius. For the Galactic 
bulge and disk components we use formulae: 




Mb 


%(jc,y,z) 


-Mb 

^Jx^+y^+z^ + bl 


<f>i/(r,z) = 

= ‘i’rf(^,y,z) = 


Md 

-Md 


4 ^ 

d 

-c 

dz 


where 




■^'^h(^,y,z} = 


-^<t>b{x,y,z) = 


dx 


d 

d~y 


^d(.x,y,z) = 


.y,z) 

-At,, 

{x,y,z)- 

-^<5„(x,y,z) 

.y,z) 


{x,y,z)- 

-|;<5„(x,y,z) 

,y,z) 


{x,y,z)- 

-§-Mx,y,z) 


xMb 


xMb 

(x2 

+ y2 +z2 


{R^ + blY 


yMh 


yMb 

(x2 

+ y2 + z2 

+ blY 

[R^+blY 


zMb 


zMb 

(x2 

+y2 +z2 

+ blY 




xMd 



:2+y2-)- 

(ad + ^Jz^ + Fd) ) 



yMd 



'■+y'^+{ad + yJz^ + F^ 


and for the Galactic halo we have: 













































On the accuracy of close stellar approaches determination 5 



Figure 2. Distribution of clones of star HIP 19946 in the original Galactic heliocentric frame. In this and all next similar plots black points are the position 
of clones in the moment of the closest approach obtained from numerical integration. Green points are the position of clones calculated with a straight line 
motion approximation. Red dot is the position of the nominal solution from numerical integration, yellow point is the proximity position for nominal star in a 
straight line motion model. Blue circle is the boundary of the Oort cloud at 0.5 pc. 


d 

dz 


'^d{x,y,z) 


zMd [aj + 

+y2j^(^ad + ^Jz^+ blj ^ 


d 

dx 


<t^h(x,y,z) 


_ ^ ^ 

Uf, f +y^ + Z^ (aj, + f +y^ + Z^'^ 

^ xMh 

af/R (cik -\-R) 


d_ 

dy 


xt>h(x,y,z) 


_ ^ ^ 

Uf, y/x^ +y^ + Z^ (aj, + 

^ yMh 

aijR{ai,+R) 


d_ 

dz 


't>h{x,y,z) 


_ zMh _^ 

ah fx^ +z~ (^ah + y/x^ +y- + z^'^ 

^ zMh 
ahR{ah+R) 


Numerical pa rameters for this model are taken from 
lirrgang et alj ( l2013h and are summarised in Table [T] 

Since we always numerically integrate these equations sim¬ 
ultaneously for a star and the Sun (to present final, heliocentric 
results), we also account for the gravitational attraction between 
these two bodies, as they can approach arbitrarily close. Our test 
calculation shows that the greatest influence of taking into account 
Sun-star gravitational interaction is for star HIP 26744. This dif¬ 
ference in the proximity distance is 0.00049 pc, less then 0.03 per 
cent of nominal value for this star. On this basis, we can state that 
interactions between Sun and stars can generally be neglected in 
such investigations. 


2.4 Principal components analysis 

The result of our calculation is always in the form of a swarm of 
10000 clones stopped at the closest proximity from the Sun. We 
have inspected all these swarms of clones in 3D and found, that they 
are typically very flat and can be conveniently presented in 2D plots 
after applying necessary rotations. To this aim we determine the 
plane of the maximum scatter using p r incipa l components analysis 
(PCA). As can be readed in Ijolliffa ( l2002h principal component 


analysis is one of the oldest techniques of multivariate analysis. 
Using this technique, we can find the largest scatter plane for our 
clones which is what we would like to present in our plots. Raw 
results of our calculations consists of Galactocentric coordinates of 
each star and the Sun. First we must go back to heliocentric frame 
so we need to subtract the calculated position of the Sun from that 
of the star. 

Then we need to construct a covariance matrix: 


cov(x,y) cov(x,z)\ 
cov(y,y) cov(y,z) (9) 

cov(z,y) cov{z,z)J 

where for example: 


( cov(x,x) 
cov(y,x) 
cov(y,x) 


t ^ I."=]ixi-x){yi-y) 

cov(x,y) =-;-= cov(y,x) 

)j — 1 

and Xi,yi,Zi - are coordinates of ith clone; x,y,z - are respective 
means of coordinates of all clones; u - is the number of clones (here 
typically 10000). 

Next we use Jacobi transformation to find eigenvalues and cor¬ 
responding eigenvectors. Then we create matrix 3x3 where in first 
column we place eigenvector corresponding to the largest eigen¬ 
value and in next columns eigenvectors corresponding to the re¬ 
maining eigenvalues. When we have this matrix with eigenvectors 
V, and vector with coordinates (x, y, z) of clone in Galactic helio¬ 
centric frame we can calculate new coordinates (x', y' and z') for all 
clones from equation [Tol 


= (Vl V2 



( 10 ) 


When we look at the results in the original Galactic helio¬ 
centric frame (Fig. [2j we see a comparable level of scattering in 
all planes, but if we express our results with respect to the frame 
obtained with PGA (Fig. O we see large difference in scatter level 
between x'y' plane and the other two. This clearly reflects the fact 
that the swarm of clones stopped at the closest proximity with the 
Sun is almost flat. 

The final result from the PGA (for plotting purposes) we ob¬ 
tain when we reduce our data to two dimensions. To reduce dimen¬ 
sions we use only two eigenvectors corresponding to the two largest 
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Figure 3. Distribution of the same swarm of clones of star HIP 19946 as in Fig.[2|but in jdy'z! heliocentric frame obtained from the xyz frame after rotations 
determined using PCA. 
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Figure 4. Distribution of clones HIP 19946 using PCA with reduction to 
2D 


eigenvalues (eq. [ID- 



We see that the final plot (Fig. |4j strictly correspond to the first 
projection in Fig.jD 

For the star used in the above examples a straight line motion 
model is a good approximation for calculating distance of closest 
approach to Sun, but in next section we show examples when this 
approximation it’s not that good. As mentioned earlier a distance 
calculated with a rectilinear motion approximation may be signi¬ 
ficantly different from a distance obtained from a numerical integ¬ 
ration. 

3 RESULTS FOR SELECTED STARS 

3.1 Selection of stars and overall results 

In this research we restricted ourselves to the stars from F1IP2 cata¬ 
logue with known radial velocities. We are interested in the closest 
stellar passages so a very small subset of these stars were chosen 
for our calculations. The selection procedure was as follows. For 
each star in the HIP2 catalogue we checked, if its radial velo city is 
available in the XHIP catalogue [Anderson & Fran^ ( 12012h . If it 


is present, we calculated the minimal heliocentric distance of this 
star (in past or in future) using the rectilinear motion approxima¬ 
tion (£);). As a short list of candidates ( 2538 objects) we selec¬ 
ted all stars having this minimal distance £); smaller than 20 pc. 
We used such a large threshold value expecting that the exact min¬ 
imal distance, obtained from a numerical integration can 

be substantially different. In our sample we did not find any star 
which Dj > 10 and < 2 pc so the criterion D; < 10 would 
be probably sufficient. Than we integrated numerically equations 
of motions of every pair, the Sun and each star from the short list, 
using full Galactic potential as described in previous section. Our 
final list consists of 40 stars with the nominal proximity distance 
smaller than 2 pc. In Table |D we present a complete list of these 
stars with their common names, nominal proximity distances and 
epochs and quality of radial velocity value. In the last column we 
describe whether the star is our new findings or its proximity was 
noticed and calculated in some earlier papers. In three cases new 
stars were independently found by us and by B-J. 

In columns 5 and 6 of TablejDwe present results (minimal dis¬ 
tance D,„i„ in parsecs and the corresponding moment of time r,„,„ 
in Myr from the present epoch) for the nominal astrometric para¬ 
meters of each star. The aim of this paper is however to estimate the 
accuracy of these results. To achieve this we replaced each star with 
10000 of its clones, obtained as described in Section lT^ Next we 
numerically integrated motion of all clones and estimated the dis¬ 
tance of the most probable proximity point and the corresponding 
epoch with their uncertainties. 

The results for all selected stars are presented in Table |D and 
here D (minimal distance) and T (the epoch of the proximity) are 
the most probable values obtained from the respective clone distri¬ 
butions. In some cases they significantly differ from nominal results 
presented in Table |D In the case of moment of time its calculating 
procedure was quite simple: T is the mean value and its uncer¬ 
tainty AT is the half of the symmetric interval covering 90 per cent 
of individual clone values. The situation is much more complicated 
with the proximity distance estimation, since their distributions are 
frequently significantly asymmetric. To present the most inform¬ 
ative result we decided to calculate three mean coordinates of the 
most probable proximity point in space (calculating simple means 
in each coordinate, i.e. x', y' and z' ) and then to calculate its dis¬ 
tance from the Sun; 

D=\lF^+y^+7^ ( 12 ) 

To estimate its uncertainty we measure the radius of a sphere 
around the proximity point in which 90 per cent of individual clone 
proximities are placed and this value is presented in TablejDas AD. 

From the inspection of the data in TablejDseveral observations 
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Table 2. List of stars with nominal proximity distance smaller than 2 pc. Third column shows quality of radial velocity quoted from the XHIP catalogue. 
Results from numerical integration of nominal data D^m and Tm,-,, are presented in fourth and fifth column, in next two columns we show for comparison 
purpose the results from the rectilinear motion approximation D/ and T;. 


HIP2 ID 

Name 

q_RV 

Dfnin [p*^] 

T„,m [Myr] 

Di [pc] 

7] [Myr] 

Earlier publications 

1392 

HD 1317 

A 

1.70 

-0.48 

1.70 

-0.48 

new 

3829 

Van Maanen Star 

D 

0.95 

-0.02 

0.95 

-0.02 

a, c, f 

12351 

GJ1049 

C 

1.93 

-0.62 

1.92 

-0.62 

a, c, e, f 

14473 

HD 19376 

A 

0.35 

-3.73 

2.37 

-3.78 

new 

14754 

GJ 127. lA 

D 

1.62 

-0.29 

1.63 

-0.29 

a, b, e, f 

19946 

BD+03 580 

A 

1.80 

-0.54 

1.81 

-0.54 

new 

21539 

CD-33 1835 

D 

1.92 

-0.14 

1.92 

-0.14 

new 

23415 

HD 32111 

A 

1.71 

5.20 

3.09 

5.16 

new 

25001 

HD 34790 

D 

1.93 

4.42 

1.88 

4.39 

a, c 

25240 

HD 35317 

A 

1.66 

-0.99 

1.63 

-0.99 

a, c, f 

26335 

HD 245409 

A 

1.54 

-0.49 

1.53 

-0.49 

a, b, c, d, e, f 

26624 

HD 37594 

A 

1.98 

-1.87 

1.95 

-1.87 

a, c, e, f 

26744 

HD 37374 

A 

1.78 

13.99 

2.26 

13.19 

a, c 

27288 

HD 38678 

A 

1.31 

-0.85 

1.30 

-0.85 

a, b, c, e, f 

30067 

HD 43947 

A 

1.78 

-0.66 

1.78 

-0.66 

a, c, e, f 

30344 

HD 44821 

A 

1.10 

-1.56 

1.11 

-1.56 

a, b, e, f 

33369 

V*BG Mon 

B 

0.85 

-5.41 

3.68 

-5.38 

new 

38228 

HD 63433 

A 

1.95 

1.35 

1.93 

1.35 

a, b, c, e, f 

38965 

HD 66589 

A 

1.79 

-1.09 

1.80 

-1.09 

new, f 

42525 

BD-t41 1865 

A 

0.82 

-0.24 

0.82 

-0.24 

new, f 

47425 

GJ 358 

C 

1.87 

-0.06 

1.87 

-0.06 

a, c, f 

54035 

HD 95735 

A 

1.44 

0.02 

1.44 

0.02 

a, b, c, d, e, f 

57544 

GJ445 

A 

1.06 

0.05 

1.06 

0.05 

a, b, c, d, e, f 

57548 

GJ447 

A 

1.92 

0.07 

1.92 

0.07 

a, b, c, d, e, f 

63721 

HD 113447 

A 

0.12 

0.13 

0.12 

0.13 

new, f 

70890 

Proxima Centauri 

B 

0.94 

0.03 

0.94 

0.03 

a, b, c, d, e, f 

71681 

CL Centauri B 

A 

0.96 

0.03 

0.96 

0.03 

a, b, c, d, e, f 

71683 

CL Centauri A 

A 

0.98 

0.03 

0.98 

0.03 

a, b, c, e, f 

75311 

BD-02 3986 

A 

1.62 

4.84 

3.55 

4.92 

a, c 

77910 

HD 142500 

A 

1.82 

3.05 

1.04 

3.07 

a, c 

84263 

HD 155117 

A 

1.11 

-6.53 

2.44 

-6.52 

new 

85661 

HD 158576 

A 

0.52 

1.88 

0.39 

1.88 

a, c, f 

87052 

HD 161959 

A 

1.91 

5.80 

6.92 

5.88 

new 

87937 

Bernard Star 

A 

1.15 

0.01 

1.15 

0.01 

a, b, c, d, e, f 

89825 

GJ710 

A 

0.29 

1.39 

0.30 

1.39 

a, b, c, d, e, f 

90112 

HD 168769 

A 

0.89 

-1.86 

0.97 

-1.86 

a, c, f 

92403 

GJ729 

A 

1.98 

0.15 

1.98 

0.15 

a, b, c, d, e, f 

94512 

HD 179939 

A 

1.99 

3.63 

2.04 

3.63 

a, c, f 

103738 

HD 199951 

A 

1.18 

-3.86 

2.36 

-3.90 

c, f 

110893 

HD 239960 

A 

1.92 

0.09 

1.92 

0.09 

a, b, c, e, f 


-Ijimenez-Torres et alJ ^201 1 

1. b -iBobvlevI ilOlOah. c -iGarci'a-Sanchez et al.l 

200 ih 

-iDvbczvnski & Kankiewicz 

Jl999h. e -iDvbczvnskil J200fl). f -iBailer-JonesI 

2014^ 


can be made. First, there are two stars, FlIP 33369 and FlIP 75311 
for which our results are completely unreliable. First star is at the 
exceptionally large distance at the present epoch (424 pc) and any 
uncertainties in its position and/or velocity significantly amplifies 
when numerically integrating its motion back to the solar proxim¬ 
ity. While its nominal proximity distance is only 0.85 pc we, basing 
on the contemporary astrometric data for this star, cannot say any¬ 
thing reliable about its real minimal distance from the Sun 5 Myr 
ago. 

For the second star, HIP 75311, the main source of the unac¬ 
ceptable uncertainty in the proximity distance is the large formal 
error of its proper motion. Both components of the proper mo¬ 
tion of this star have uncertainties on the level of 300 per cent in 
the HIP2 catalogue. The nominal proximity distance is 1.62 pc, 
the most probable value equals 5.71 and both should be treated as 


highly unreliable. Both these stars are marked with ’?’ in Table[3as 
well as five more stars due to their large proximity position errors. 
There is also one additional star marked in the same way despite 
of its small error, namely HIP 63721. The reason for this is that 
we are aware of its highly unreliable parallax in HIP2. Instead of 
;r = 217 mas included in HIP2 othe r sources present parallax valu e 
well beyond 5 mas (see for example [Fabricius & Makar^ ( |200(]|) ). 
HIP2 value is also discarded in XHIP catalogue. 

There are also 12 stars (marked with just after their HIP 
numbers in Table O in our final list which are now close to the 
Sun and simultaneously close to their proximity epoch. For these 
stars a proximity uncertainty is practically equal to their current as¬ 
trometric position and velocity uncertainties and therefore small or 
very small. In our calculations we integrated them only over ex¬ 
tremely short time intervals and we did not observe any uncertainty 
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Table 3. Estimated minimum distance R from the Sun and corresponding moment of time T for 40 stars that nominally can came closer than 2 pc. A dash 
marks cases where we cannot obtain reasonable values. 


HIP ID 

D[pc] 

AZ)[pc] 

T [Myr] 

AT [Myr] 

HIP ID 

D[pc] 

AD [pc] 

T [Myr] 

AT [Myr] 

1392 

1.70 

0.16 

-0.48 

0.02 

47425 * 

1.90 

0.45 

-0.06 

0.01 

3829 * 

0.95 

0.05 

-0.02 

<0.01 

54035 * 

1.44 

0.01 

0.02 

<0.01 

12351 

2.26 

1.47 

-0.71 

0.45 

57544 * 

1.06 

0.03 

0.05 

<0.01 

14473 ? 

0.22 

7.84 

-3.78 

0.74 

57548 * 

1.92 

0.04 

0.07 

<0.01 

14754 

1.64 

0.27 

-0.29 

0.04 

63721 ? 

0.13 

0.04 

0.14 

0.03 

19946 

1.91 

1.32 

-0.55 

0.13 

70890 * 

0.94 

0.03 

0.03 

<0.01 

21539 * 

1.97 

0.49 

-0.14 

0.03 

71681 * 

0.96 

0.09 

0.03 

<0.01 

23415 ? 

1.65 

7.82 

5.35 

1.43 

71683 * 

0.98 

0.04 

0.03 

<0.01 

25001 

2.07 

1.54 

4.67 

1.70 

75311 ? 

5.71 

- 

5.23 

- 

25240 

1.66 

0.30 

-0.99 

0.05 

77910 

1.94 

1.59 

3.15 

0.90 

26335 

1.54 

0.09 

-0.49 

0.01 

84263 ? 

1.20 

4.64 

-6.59 

1.07 

26624 

1.98 

0.25 

-1.87 

0.08 

85661 

0.52 

0.62 

1.88 

0.11 

26744 ? 

1.71 

4.21 

14.12 

2.30 

87052 ? 

1.81 

5.87 

5.90 

1.27 

27288 

1.31 

0.10 

-0.85 

0.06 

87937 * 

1.15 

0.01 

0.01 

<0.01 

30067 

1.78 

0.11 

-0.66 

0.01 

89825 

0.29 

0.35 

1.39 

0.08 

30344 

1.11 

0.23 

-1.57 

0.11 

90112 

0.91 

1.31 

-1.87 

0.20 

33369 ? 

51.00 

- 

-4.96 

- 

92403 * 

1.98 

0.04 

0.15 

<0.01 

38228 

1.95 

0.16 

1.35 

0.03 

94512 

2.06 

2.33 

3.65 

0.46 

38965 

1.87 

1.16 

-1.10 

0.23 

103738 

1.18 

0.69 

-3.87 

0.30 

42525 

0.99 

0.96 

-0.25 

0.09 

110893 * 

1.92 

0.04 

0.09 

<0.01 
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1 



-1 


-2 -'-1-^-1-'-1-'- 

-101 23 

X' [pc] 

Figure 5. HIP 25240 - the example of a compact swarm of clones resulting 
from a cun'ent distance of this star which is relatively close to the minimum 
value (nominally 1 Myr ago). 


amplification. The example of such a result is presented in Fig|5] 
Any further discussion of the accuracy of our results for these stars 
seems unnecessary. 


3.2 Examples 

The remaining 20 stars present a wide variety of proximity dis¬ 
tance accuracy estimations. We will discuss several representative 
examples in following subsections. 


3.2.1 HIP 89825 ( Gliese 710) 


This star 
neighbourhood, 


a well known future 
see for example 


visitor in the solar 


Miillari & Orlov 


Dybczynski & Kankiewic j ( Il999l) : iGarcia-Sanchez et al.l (2001) 
Dvbczvhskil (l200fil) . According to the latest estimation of its 


1996); 


astrometric paramet ers presented in H 1P2 catalogue and using 
Vr = —13.8 kms^* dGontcharovI ( l2006h . used also in XHIP) we 
obtained the proximity distance (both nominal and mean) as small 
as 0.29 pc. As it is clearly depicted in Fig|^ the straight linear 
motion approximation works very well in this case. The swarm 
of clones of this star is however significantly dispersed, mainly 
due to large uncertainties of proper motions of Gliese 710. If 
one had used an older value of its radial velocity: -23 kms^* 
dWilsonll 19530. later con firmed for example in the Palomar/MSU 
survey dPeid et al.|ll995[) , the minimal distance would reduce even 
further, down to 0.17 pc and almost all clones of this star would be 
situated closer than 0.5 pc from the Sun i.e., the widely accepted 
radius of the outer Oort cloud. 


3.2.2 HIP 14473 (HD 19376) 

This star is our new finding and the only one in our list for which 
an estimated proximity distance in the past (3.78 Myr ago) is sig¬ 
nificantly smaller than the adopted radius of the Oort cloud. We 
estimated this minimal distance to be D =0.22 pc with the prox¬ 
imity position uncertainty AD =7.84 pc (see Fig|7ll. The reason of 
such a large error of our estimation is a large uncertainty of proper 
motion of this star. If one obtain much more accurate value (for 
example from Gaia mission) this star should be considered as a ser¬ 
ious candidate for the stellar perturber of cometary motion in the 
past. 


3.2.3 HIP 103738 (HD 199951) 

Here we have a good example of the significant difference between 
the exact numerical integration in the Galactic field and the rectilin- 
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Figure 8. Distribution of clones of HIP 103738 
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Figure 7. Distribution of clones HIP 14473 
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Figure 9. Distribution of clones HIP 87052 




ear approximation (see Fig[^. Using the numerical integration we 
obtained the minimal distance smaller by a factor of two. This star 
is now relatively distant ( about 70 pc) but all its positional and kin¬ 
ematic parameters are known with such a good accuracy that even 
long term numerical integration (for almost 4 Myr) do not down¬ 
grade the precision of the result significantly. 

3.2.4 HIP 87052 (HD 161959) 

This star is another one among 10 stars recognised in this paper for 
the first time as possibly visiting the solar vicinity. The reason for 
this novelty is that the straight line motion approximation (green 
points in Fig. 0 gives the proximity distance for this star as large 
as almost 7 pc! The distance to the most probable proximity point 
from numerical integration obtained by us is only 1.81 pc but its 
dispersion is large so we estimated rather large error for this value, 
namely ±5.87 pc. This is the consequence of a large uncertainty of 
proper motion for this star. This is also present when we compare 


nominal proximity points obtained on the basis of different astro- 
metric catalogues, as it is depicted in Fig.llOl 

3.2.5 HIP 84263 (HD 155117) 

In this case the minimal distance obtained from a numerical integ¬ 
ration reduces to 1.20±4.64 pc comparing with 2.44 from a rec¬ 
tilinear model. This star has very small proper motion with rather 
big relative errors which together with the long numerically integ¬ 
rated time interval (over 6 Myr) gives a widely spread swarm of 
clones, with significant number located inside the assumed Oort 
cloud sphere (see Fig. [ID- 

3.2.6 HIP 77910 (HD 142500) 

Using numerical integration with full Galactic potential we almost 
double the minimal distance of HIP 77910 from the Sun, obtain¬ 
ing D = 1.94 ± 1.59 pc (black points in Fig. d. This makes its 
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Figure 10. Proximity point positions of HIP 87052 obtained from differ¬ 
ent catalogue data. This plot is a projection on the plane peipendicular to 
present heliocentric direction to this star. Different symbols correspond to 
different source catalogues and different colours correspond to three differ¬ 
ent radial velocity values used in this plot. 
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Figure 11. Distribution of clones of HIP 84263 


gravitational influence on Oort cloud comets weaker but since its 
estimated mass is 2 Mq it should still be considered as a potential 
stellar perturber. 


3.2.7 HIP 57791 (HD 102928) 

HIP 57791 is an extra example (outside our star list but for ex¬ 
ample included bv lGarcia-Sanchez et al.|[l99^ of completely dif¬ 
ferent results from the rectilinear approximation and numerical in¬ 
tegration. In this case a full model moves out the proximity point 
from 1.8 pc to 3.7 pc from the Sun (see Fig. I13t what makes this 
star rather improbable stellar perturber for the Oort cloud comets. 
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Figure 12. Distribution of clones of HIP 77910 
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Figure 13. Distribution of clones HIP 57791 

4 ’LAST MINUTE’ DISCUSSION 

Just in the moment we had our pa per ready for subm ission we 
learned about a very recent paper bv lBailer-JonesI ( I20l4) . He used 
almost identical data and methods and obtained results similar to 
ours in many respects. However there are also several differences 
in our approach and as a consequence his results are in many cases 
significantly different so it is worth to study in detail each indi¬ 
vidual star case. 

In his Table 3 B-J describes 65 close (closer than 2 pc from 
the Sun) stellar passages. However, for many stars he included sev¬ 
eral variants (based on different stellar data) so in fact he presented 
stellar proximities for 42 individual stars. This is very close to 40 
stars in our Tables [2] and [3] but surprisingly we have only 27 stars 
in common. There are two main reasons for this discrepancy. First, 
B-J presents his results in significantly different manner. He calcu¬ 
late simple mean of all proximity distances for all clones with 90 
per cent confidence intervals and in many cases this method some¬ 
what overestimates the distance we should expect. The heliocentric 












On the accuracy of close stellar approaches determination 11 


distance is obviously positive so its distribution is often signific¬ 
antly asymmetric. Moreover mean value is rather sensitive to out¬ 
liers and due to the limitation to positive values almost all outliers 
are placed at one side of the mean. For this reason we decided to 
present the distance to the most probable proximity position, as it 
is described in detail in Sec l3.ll Since B-J selected stars for his 
published Table 3 using only mean distances several objects were 
shifted out to his Internet large table. The most spectacular example 
is the star HIP 14473. As it is shown in Fig|7] spatial positions of 
proximity points of individual clones are widely spread. B-J cal¬ 
culated the mean of all this distances and report the result (in his 
extended, Internet table) to be 6 pc. In our approach we report the 
proximity distance as 0.22±7.84 pc (of course the second number 
is the uncertainty in proximity point position, not the formal error 
in the proximity distance which must be positive). 

The second reason for several discrepancies between B-J and 
our results is his usage of few questionable radial velocities from 
RAVE catalogue. The great majority of velocities taken by B-J from 
RAVE catalogue seems to be obviously erroneous, in some cases 
greater than the escape velocity from our Galaxy. B-J also com¬ 
mented all this cases as uncertain. In our investigation we did not 
use this velocities. 

There are also several other differences in the methods (e.g. 
the Galaxy potential model) and result presentation, including the 
main purpose of the work: we aimed at studying and presenting 
the accuracy of proximity calculations while B-J concentrate on 
completing the list of closest stellar approaches. He also discuss in 
detail many problematic cases which helped us to reveal the source 
of some discrepancies, e.g for HIP 85605. 

As it concerns the influence of the Galaxy pote ntial model, 
B-J used the Dauphole and Colin model i ll9951 : 1 19961) and states: 
“While better models may now exist, the results are not very sens¬ 
itive to the exact choice”. The influence of different Galaxy models 
for close stellar approaches to the Sun was recently discussed by 
Ilimenez-Torres et aP ) l201ll) and their conclusion is that for more 
distant stars the results are sensitive to the choice of the model. 
They present a large table comparing results for different models 
(their tab.3) but erroneously instead of presenting t heir o wn res¬ 
ults they copied results form iGarcfa-Sanchez et al. ) l200ll) . Since 
we used a very recent model by Irrgang et all |20H ) we are inter¬ 
ested in such a comparison. 

B-J kindly made available for us his nominal results (Bailer- 
Jones, 2015, personal communication). Thanks to that we were able 
to study the difference between stellar results for two models in use 
(including slightly different initial Sun position and velocity and 
marginal difference in the Galactic frame orientation). We found 
moderate differences in the nominal closest stellar approach dis¬ 
tances for several objects. In Table|4]we present selected examples 
of such differences (for the purpose of this comparison we used 
exactly the same stellar initial data as B-J). 

In the following sections we discuss in detail all discrepancies 
we found in the stellar close approaches lists of B-J and ours. 


4.1 14 stars present in B-J but omitted here 

There are two groups of stars from B-J list missing in our paper. 
First is the result of using RAVE radial velocities (not used by us) 
and the second comes from the utilisation of astrometric parameters 
(mainly proper motions) from the XHIP catalogue, where HIP2 as¬ 
trometry was replaced by some other, frequently taken from HIPl 
or Tycho-2. 

From the first group (10 stars), in our opinion, only star 


Table 4. The examples of the minimal star-Sun distances (nominal) cal¬ 
culated with different Galaxy potential model. is calculated by B-J 
(Bailer-Jones, 2015, personal communication) with the Dauphole and Colin 
model h995l:ll99el) Pj is calculated by us with the model used in this 
paper lirrgang et alj|20l3) . Object names consist of the HIP number aug¬ 
mented with a letter indicating the source of a radial velocity, detailed ex¬ 
planation can be found in B-J. 


object 


Dlin [pc] 

14473g 

0.072 

0.375 

14473p 

0.095 

0.353 

75159r 

0.347 

0.251 

23415p 

1.486 

1.715 

84263x 

1.419 

0.802 

84263p 

0.723 

1.112 

84263g 

0.731 

1.117 

87052p 

0.394 

1.914 

87052x 

1.409 

2.136 

91012r 

0.446 

0.350 


HIP 87784 can be treated as having close approach to the Sun (as¬ 
suming that the RAVE radial velocity of -66.30±3.60 km s^* found 
by B-J will be confirmed). Also HIP 75159 might be considered 
as real (but highly uncertain) solar neighbourhood visitor. The 
rest (HIP 23311, HIP 4I3I2, HIP 53911, HIP 55606, HIP 9I0I2, 
HIP 100280, HIP 104256 and HIP 104644) are results of using un¬ 
acceptable radial velocities. 

Coming to the second group: 

• HIP 34617 is the double system and proper motions in HIP-2 
are not necessarily the best, so probably B-J result, using Tycho-2 
values, is better. 

• HIP 85605 is also a double system and erroneous paral- 
lax for this star exi sts in Hipparcos, as pointed out by SIMBAD 
dWenger et al.ll200^ database comment (“Parallax and proper mo¬ 
tion are not compatible with CCDM J17296-1-2439A; the large 
proper motion and parallax of this star in Hipparcos (;r=202mas, 
jt=362mas/ vr) is most likely an artefact” ). The similar opinion is 
expressed in ipabricius & MakarovI ( l2000l) . 

• HIP 86961 and HIP 86963 are members of the triple system 
and according to lHenrv et al.l ( l2006l) both these stars have signific¬ 
antly erroneous parallaxes in Hipparcos catalogues. 


4.2 12 stars present in this paper but omitted in B-J 

Apart from some differences in Galaxy potential model described 
above the main reason for some our stars missing in B-J list is 
presenting the mean value of the closest Sun-star distance. In many 
cases this results in larger values than presented by us, which 
includes stars HIP 1392, HIP 14473, HIP 19946, HIP 23415, 
HIP 25001, HIP 26744, HIP 77910, HIP 84263 and HIP 87052. 
For all these 9 stars the nominal proximity distance calculated by 
B-J is below the 2 pc threshold but his mean distance is larger and 
he excluded them from his published list. 

There are three more stars in our list that are omitted by B-J. 
First is star HIP 21539 for which nominal proximity distance for 
data taken from XHIP catalogue is 1.92 pc but uncertainty of its 
radial velocity is undetermined (see Sec. l2.2t and probably for that 
reason B-J removed this star from his sample. Next two stars are 
HIP 33369, our nominal result equals 0.85, and HIP 75311 - HIP2 
astrometry combined with the radial velocity from Pulkovo cata¬ 
logue gives nominal proximity distance of 1.62 pc. These two stars 
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have enormously large uncertainty of the proximity point position 
so their mean distance values are larger then 10 pc (and therefore 
omitted by B-J), but their nominal solutions are below the 2 pc 
limit. 


4.3 Other difficult cases 

As we have all nominal solutions kindly made available by B-J 
(also those corresponding to his large, Internet table) we carefully 
compared his results with ours. Few problematic cases we found 
for stars with the nominal proximity distance smaller than 2 pc: 

• For HIP 1647 B-J found a radial velocity of -86.8 ± 28.8 
kms^* in RAVE catalogue and obtained the closest nominal dis¬ 
tance of 1.979 pc. While in our model this distance is even smaller, 
1.881 pc we think this result should be discarded. The reason is 
that HIP 1647 have a well established (since 1928) radial velocity 
of - 1-12 kms^*, i.e. in opposite direction. 

• For HIP 10332 B-J have obtained a minimal, nominal distance 
of 1.793 pc basing on XHIP data. While XHIP authors treat this 
star astrometry as problematic and incorporate (a little bit inconsist¬ 
ently) its parallax from HIPl and its pro per motion from T ycho-2, 
the problem seems to be more serious. In lLuck et al.l ( 1201 ll) the au¬ 
thors obtained the heliocentric distance for this star to be over 4 kpc. 
HIP 10332 (V*UX Per) is a classical 5 Cepheid type star with the 
period over one year. There was well known problem in Hipparcos 
mission with obtaining correct parallaxes of classical cepheides and 
this star is most probably one of the problematic Hipparcos results. 

• B-I used also XHIP data for HIP 24670 obtaining the nominal 
proximity of 1.729 pc. Proper motions of this star are very small 
with large uncertainties in all catalogues, which connected with its 
large current distance (over 160 pc according to HIP2) makes the 
uncertainty of the proximity distance very large. Mean value from 
B-J equals 5.3 pc. 

• HIP 32475 according to B-J has the close approach nominal 
distance of 1.725 pc basing on XHIP data, where TYCHO-2 proper 
motions are included. This star is in fact a close binary (CCDM 
J06467-I-0822AB), unresolved in Hipparcos and Tycho catalogues. 
If we treat TYCHO-2 proper motions as superior to HIP2 this B-J 
result is fully valid. 

• For HIP 75807 B-J found a radial velocity of -37.40±1.20 in 
RAVE catalogue which, combined with HIP2 astrometry gives a 
nominal proximity distance of 1.461 pc. The problem is that the 
HIP2 astrometry for this star (especially its parallax) is evidently in 
error. In the same RAVE catalogue one can find the spectrophoto- 
metric parallax for this star equal to 2.55 ma s, which is highly con¬ 
sistent with HIPl result, further corrected bv ipabricius & Makarov! 

( l200(]h . 


5 CONCLUSIONS 

The main purpose of this investigation was to study the accuracy 
of determination of the proximity distance and epoch for stars that 
can come ( in past or in future) close to the Sun and act as stellar 
perturbers of the long period comet motion. Our study is based on 
the best available da ta in the pre -Gaia era, mainly from the HIP2 
dvan Leeuwetill201 ll) and XHIP d Anderson & Francisll2012h cata¬ 
logues. We used a numerical integration in the rectilinear, Galacto- 
centric coordinates taking into account th e full gravitational p oten- 
tial of the Galaxy in its modern shape dirrgang et alj|201^ . For 
the main purpose (the accuracy assessment) we utilised covariance 


matrices of astrometric data which is included in the HIP2 cata¬ 
logue. Each star was substituted by a swarm of 10 000 of its clones 
according to the respective covariance matrix and radial velocity 
dispersion. Than we used Principal Component Analysis to present 
the distribution of the swarm at the proximity epoch. 

Limiting our self to the nominal proximity distance of 2 pc 
we formulated a list of 40 stellar close visitors and among them 
10 stars were our independent new findings (three of them were 
independently pointed out by B-J). We analysed each case in detail 
and concluded that for more than 50 per cent of stars in our list the 
accuracy is good or very good because these stars are now close to 
their proximity epoch. We showed several examples of moderate 
accuracy for more distant (in space and in time) stars. 

We also showed that the linear approximation method can lead 
to large errors in many cases, especially for more distant stars, 
where the curvature in their motion induced by a Galaxy gravity 
cannot be ignored. 

We also concluded that for two stars (HIP 33369 and 
HIP 75311) the extremely large formal errors makes our proximity 
distance results completely unreliable. We suspect that the result 
for HIP 63721 is also erroneous due to its probably false parallax 
in the HIP2 catalogue. We also marked next five stars ( HIP 14473, 
HIP 23415, HIP 26744, HIP 84263 and HIP 87052) as having un¬ 
acceptably large estimated errors. The first of them might be an im¬ 
portant perturber of the long period comets motion since its nom¬ 
inal proximity distance was as small as 0.22 pc about 3.78 Myr ago. 
The problem is that the proximity position uncertainty (AZI = 7.84) 
is unacceptably large due to the astrometry uncertainties. 

We are convinced that the incoming advance in stellar astro¬ 
metry caused by Gaia mission can significantly improve the accur¬ 
acy of nearby stars passages close to the Sun. The most important 
part is the precision of proper motions and parallaxes, augmented 
with (mainly ground based) precise radial velocity measurements. 
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